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We present results of a quenched calculation of the heavy-light and quarkonium spectrum using the tadpole 
improved clover action. We resolve completely the triplet \ P-states in quarkonium systems, and obtain evidence 
for fine structure of the heavy-light P-states. Approximate scaling of the hyperfine splittings is observed, producing 
results that are significantly below experiment. 



1. Introduction 

Charm physics continues to pose difficulties for 
lattice QCD simulations. The systems are signif- 
icantly relativistic causing problems for the 
NRQCD approach, and have arriQ ~ 0(1) caus- 
ing significant discretisation effects in the heavy 
Wilson quark approach. We present the results 
of a simulation using the tadpole improved clover 
action and perform the analysis using the Fermi- 
lab H interpretation of the heavy Wilson quark 
approach. 

2. Simulation Details 

The simulation was performed using 499 
quenched gauge configurations on a 16 3 x 48 lat- 
tice. Five heavy quark masses and three light 
quark masses, detailed in Table [l], were simu- 
lated using the tadpole improved clover action, 
with t*o = 0.8778 from the average plaquette, 
and C S w = 1.47852. Both local and fuzzed [Q 

Table 1 



(3 = 6.0 simulated kappas 



K 


aMps 


Fuzzing Radius 


0.13856 


0.228(2) 


6 


0.13810 


0.293(1) 


6 


0.13700 


0.4135(10) 


6 


0.13000 


0.9283(7) 


3 


0.12600 


1.1618(6) 


3 


0.12200 


1.3755(6) 


3 


0.11800 


1.5751(6) 


3 


0.11400 


1.7644(6) 


3 



operators were generated at source and at sink. 
(Local) covariant derivative sources in each of the 



spatial directions were used for the k = 0.12600 
quark corresponding to K c h arm , allowing the op- 
erator for a 3 P2 state to be created for combina- 
tions involving k = 0.12600 with each of the other 
masses. The operators used are given in Table |. 

Table 2 



Meson operators 



State 


J PU 




Operators 




0"+ 








1~ 




irnip 


x Pi 


1+- 




lp<Hjip 


3 Po 


0++ 




ipip 


3 P1 


1++ 






3 P 2 


2++ 




— 7j Aj}ip E rep 








+ T rep 



An extensive analysis of correlated double and 
single exponential fits to various smearing combi- 
nations was carried out, and the optimal fitting 
approach selected for each channel in the light- 
light, heavy-light and heavy- heavy sectors inde- 
pendently. 

3. Fine Structure 

We obtain a signal for the fine structure of the \ 
triplet of P-states in quarkonium, as illustrated in 
the effective mass plots in Figure [|. In the heavy- 
light case we obtain a signal for the splitting be- 
tween the 1 + and + correlators, Figure |, fitted 
using a single exponential model for the fuzzed- 
fuzzed combination. The consistency of both 
double and triple exponential multi-correlator fits 
was checked. Cross correlations between the 1 ++ 
and l" 1 operators showed that there was signif- 
icant mixing, as expected since the jn g ht basis is 



2 




Figure 1. Quarkonium \ triplet for the (k — 0.12600, k = 0.12600) combination. Double exponential fits were 
performed to the fuzzed-local and local-local correlators simultaneously for the ++ and 1 ++ states, while a single 
exponential fit was performed to the derivative based operator for the 2 ++ state. The three states are clearly 
resolved. 



Figure 2. Heavy- light fine structure for the (kq — 
0. 12600, n q = 0.13810) combination; fit is to timeslices 
7-13 of the fuzzed-fuzzed correlator. 

— , , , | , , | r- 

+ 



the physical basis in the ttiq — > oo limit. However 
we could not resolve the two physical 1 + states, 
and treat the ground state for each of the 1 + cor- 
relators as the lower lying 1 + state, assumed to 
have jiight = \. 



4. Spectrum at (3 = 6.0 

We use the dispersive mass of pseudoscalar 
mesons defined by E(p 2 ) — mi + + Cp 4 , as 
the definition of the heavy quark mass in extrapo- 
lations of the spectrum. We found that linear ex- 
trapolations in the inverse heavy quark mass gave 
a very acceptable x 2 /dof . For the heavy-light sys- 
tems, linear chiral extrapolations were performed 
to K cr ;t and to ^strange, followed by linear extrapo- 
lations in the inverse heavy-strange pseudoscalar 
mass to the D „ and _B, masses. 



For the 3 i-2 state, where only combinations 
of propagators involving the k — 0.12600 quark 
(with another) could be formed from the available 
data, the non-degenerate combinations were used 
to extrapolate to the physical meson masses, in- 
troducing a small correction in extrapolations to 
the J/ip system. However, the extrapolations of 
the 3 P2 splittings to T are not well under control. 
The results obtained using the string tension, m p 
and the quarkonium S — P splitting to set the 
scale, and the kaon mass to fix K s trange, are tabu- 
lated in Table |Q and Table |. 

5. Scaling Behaviour 

Comparison to a calculation using 220 configu- 
rations at — 6.2 on a 24 3 x 48 lattice [|| with the 
same action, allows some estimate of the scaling 
behaviour to be made. We plot the lattice spacing 
dependence of the charmonium and D s hyperfine 
splittings in Figure [| 

Near-scaling behaviour is seen with both the 
string tension and with m p used to set the scale. 
Scaling is not seen with the quarkonium S-P split- 
ting. However, the lack of scaling is only a ler 
effect. NRQCD calculations [Q have found the 
T S-P splitting scaling well with m p , at about 

1 There is a systematic uncertainty of order 50 MeV 
in the values for the heavy-light S-P splitting since the 
experimental values are for the Jiight = f doublet, while 
we calculate Jiight = \ states. For the heavy-light P-states 
we have adopted the nomenclature used for the Jiight = \ 
doublet in the Kaon system by the PDG. 
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Table 3 

(3 = 6.0 heavy-light mass splittings (MeV) 



Scale 


VK 


M p 


S-P 


Expt 


D* 


-D 


110(7) 


106(8) 


129(10) 


142 


D* s 


-D s 


99(5) 


95(6) 


115(9) 


144 


D 8 


-D 


98(5) 


96(5) 


107(6) 


105 


Di 


- D 


540(30) 


530(30) 


600(30) 


459 


Dsi 


-D s 


494(18) 


480(20) 


545(20) 


460 


Di 


-D* 


45(20) 


44(20) 


47(25) 




D sl 


~D* s0 


57(12) 


56(11) 


64(13) 




B* 


- B 


41(9) 


39(10) 


59(11) 


46 


B* s 


-B s 


40(6) 


38(7) 


57(8) 


47 


B s 


- B 


90(6) 


88(6) 


110(7) 


91 


B 1 


-B 


490(30) 


480(30) 


590(40) 


419 


B s \ 


-B s 


440(20) 


430(20) 


540(30) 


446 


Si 


-B* 


50(20) 


50(20) 


60(20) 




B sl 


-B* s0 


43(11) 


42(10) 


54(12) 





Scale 


VK 


M p 


S-P 


Expt 


J/V- 




68(2) 


65(2) 


79(4) 


117 


'Pl- 


-s 


418(13) 


408(16) 




458 


Xc2 ~ 


Xcl 


81(28) 


78(27) 


93(33) 


46 


Xcl ~ 


XcO 


51(7) 


49(7) 


59(11) 


95 


Xc2 - 


XcO 


133(28) 


128(28) 


153(35) 


141 


T- 


Vb 


22(1) 


21(1) 


31(3) 


40 




-S 


366(17) 


358(20) 




460 


Xb2 - 


Xbi 


43(30) 


42(30) 


56(32) 


21 


Xbi - 


Xbo 


26(9) 


25(9) 


33(10) 


32 


Xb2 ~ 


Xbo 


65(30) 


63(27) 


85(32) 


53 



30% below its experimental value. This suggests 
that ultimately scaling with respect to the S-P 
splitting will be seen at values of the hyperfinc 
splittings above those with m p and the string ten- 
sion. Likewise we find the D and T hyperfinc 
splittings to show approximate scaling below ex- 
periment with m p and the string tension. After 
extrapolating the heavy-light results to B and B s , 
the errors are such that our values are consistent 
with experiment. 

6. Conclusions 

We resolve completely the triplet of \ states 
in the J ftp system with a relativistic action, and 
obtain evidence for fine structure in heavy-light 



Figure 3. Scaling Behaviour of hyperfine split- 
tings. 



J/^ Hyperfine Splitting 
« M„ Scale 



, Expt 



VK Scale 
+ S-P Scale 



Table 4 

(3 = 6.0 heavy-heavy mass splittings (MeV) 



D s Hyperfine Splitting 



Expt 



« M p Scale 
o VK Scale 
+ S-P Scale 



0.05 0.10 

a (fm) 



0.05 0.10 

a (fm) 



systems. We find that both the quarkonium and 
D hyperfine splittings scale well with both the 
string tension and with m p , lying significantly be- 
low their experimental values. We find that our 
results using the quarkonium S-P splitting do not 
scale well. This may be due to the use of only 
local propagators in the quarkonium calculation 
at f3 = 6.2, (which plateau at large t, making the 
plateau identification and statistical noise toublc- 
some), or to discretisation effects. Further calcu- 
lations (better smearing at (3 — 6.2 and possibly 
at higher (3) are required to resolve which is the 
case. 
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